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ABSTRACT 


A theoretical model is derived to calculate the specific acoustic impedance of the 
absorptive material covering the walls of a cavity. This model will allow the exper- 
imental determination of the specific acoustic impedance from the measurement of the 
reverberation time in a water-filled cavity. The model assumes a wall of low absorption. 


It can not be used for rigid or pressure release walls and grazzing incidence 1s excluded. 
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I INFRODUCTION 


A. MOTIVATION 

The design of coatrng materials that reduce the reflection of water borne sound is 
an miportant problem in underwater acoustics. Since itis difficult to theoretically predict 
the complex reflection coefficient of such a material, the reflection properties must be 
determined experimentally. 

At high frequencies it 1s possible to measure the reflection coefficient of a slab of 
material by measuring the pressure after reflection of a narrow beam of sound tncident 
ata given angle. liowever, this method becomies nore difficult to use as the frequency 
of the sound 1s lowered. 

Determining the acoustic impedance of a material by measuring the standing wave 
pattern in an “impedance tube”, while the standard procedure m air borne sound, is very 
difficult for water born sound, because of the impossibility of approaching “rigid condi- 
tions” for the tube wails. This is because the “pe” for water is approximatcly 5,000 times 
higher than that for air and is much closer to that of any wall material. 

Therefore it is attractive to investigate the possibility of determining the acoustic 
unpedance of a material from measurements of the reverberation ume in a water-filled 


cavity with walls lined with the material. 


Parone CrLiVE 
To derive a theoretical model which will allow the experimental measurement of 
acoustical propertics of absorption materials from: measurements of reverberation ume 


In a water filled cavity with wall ined with material. 


C. METHOD 

The model rs based on the normal mode theory in cavities. 

Ht will be shown that the reverbcration time m a cavity can be related to the 
temporal absorpuon coefficient of an acoustic field. from this fact, a relation between 
this coefficient and the complex specific acoustic impedance (resistance and reactance) 


will then obtained, under several assumptions and limitations. 


D. REMARKS 
No references were found using the proposed approach. In addition most of the 


pubhshed articles in this field are for air as the medium in the cavity. 


Il, LITERATURE REVIEW 


This ts a summary of published articles which arc, in one way or the other, related 
to the subject. These summaries will give some idea about the accepted definitions and 


approximations relevant to this problem. 


A. CHAMBER FOR REVERBERANT ACOUSTIC POWER MEASUREMENTS IN 
AIR AND WATER [REF.]] 

This was the only paper found in the Journal of the Acoustics Society of America 
where experiments made both in air and in water were reported, (in opinion corrob- 
orated by the authors). 

The motivation for this work was the interest in obtaining sound-power measure- 
ments directly, and the idea that it might be possible with the help of a reverberant 
acoustic chamber. It is stated that if the sound field consists of many randomly excited 
nodes (diffuse field), the radiation impedance of the tank is not affected by the fluid re- 
action, because the randomly phased fluid pressures cancel. 

One of the major lessons to take out from this work 1s the fact that the reverberation 
time in the water, in a similar enclosure, is smaller than that of the air by a factor of ten. 
Yhe author attributes this to the higher mean collision rate and higher sound speed in 
water. It 1s also noticeable that the air chamber gives a more uniform sound speed. This 
uniformity 1s Important when the source radiation is concentrated in narrow frequency 
bands. 

The fluid loading factor 1s given by y = p,c,/p,h2af, where p, and c, are the density 
and the sound speed in water, p, and / the density and the thickness of the wall, and f 
the frequency. If y<1, the panels are very massive (hard walls) and specularly reflective; 
if y>1, the walls are very soft and specularly reflect sound as free surfaces. When y=] 
then the acoustic impedance of the wall is complex and the reflection is no longer 
specular, this the resistive and reactive parts of the wall impedance are function of the 


angle of incidence, this gives a more diffuse sound field. 


B. REVERBERATION TIME, ABSORPTION, AND IMPEDANCE [REF.2] 
Dowell presents a rigorous theoretical model to calculate reverberation time in a 


room in terms of the impedance of absorption materials on the wall. 


tJ 


tlis approach is based on the couphng of the normal modes of a rigid wall rooin, 
through the damping due to wall absorption. 

The sound pressure in the room 1s — by the summation of the contribution of 
ach normal mode p= p,c? a, jai Wiiere p, and ¢, are the micdiuim density and specd 
of sound, respectively, P, the coefficients in a normal mode expansion for pressure, F, 
for a rectangular room is cos(u,7x/L,) cos(m,my/L,) cos(a,z72/L,) and Af, =e,¢,e, (for 

=1 for n=0 and ¢,=1/2 otherwise). 7’, must satisfy the differential equation 
Pa. wi P + ( Pee oe: /AC,=0 with A, being the area covered with absorption 
materials, w, the natural frequency Of the mode and ©. = {(FF/2,) GAA and 2 the 
impedance of the matenial. 

From the assumption that the sound field 1s diffuse (1.c., uniform all over the room), 
the initial conditions for p. p, P, and P, can be calculated. In this case if we assume large 


impedance on the walls (z,/p,¢c,>1), all the coupling between modes are weak, and as 
Ve: 


71 


= PP =0 at 1=0 for all modes except n = 0, the calculations can be made only con- 
sidering this mode. Other approximation can be made if we consider the initial potential 
and kinetic energy on each mode as being the same, then we can find the initial cond)- 
Hens for, respectively, 7 and 7". 

In the case of low absorption walls...e., large impedance, the damping coupling co- 
efficient is small — Dowel) “ound ~— the modal réverberauen, “time 
T,=6in lor CE rs dala 


dominates, the reverberation ume rs T.. 


n 


aeC cope Or 7. wien value 1s 1 Zao est arsmivlic niode 

Dowell claims that in the standard hterature “the reverberation time 1s related to a 
random-absorption coefhicient, the random absorption coefficient is related to the 
normal-absorption coefficient, and. finally. the normal-absorption coefficient 1s related 
to impedance” [Ref. I: p. 183], it does not consider room geometry, and considers uni- 
form absorption. Iie considers his model as a generalization of the previous methods, 


since it computes reverberation time directly from wall impedance and geometric factors. 


C. SABINE REVERBERATION EQUATION AND SOUND POWER 
CALCULATIONS [REF.3] 

Various experiments show that the reverberation time changes with the amount and 
distribution of the absorption material in the cavity. This conflicts with the assumptuon 
commonly stated in the classical theorv. 

Young savs the measurement of the decav time instead of the reverberation time is 


easier and more direct. It is also advantageous when several absorption materials are 


or] 


tested on the same enclosure, because the results are obtained by comparison and the 
use of the decay time eliminates the eventual effects of any bulk absorption in the me- 
dium. 

In order to obtain a diffuse sound field, required for the measurements, one must 
wait a while after the source 1s turned on, this 1s to give time for several reflections to 


OCCUR, 


D. REVERBERATION TIME IN ENCLOSURES [REF.4] 

In this paper, the effect of the dependence of the sound absorption coefficient on the 
incident angle was tested. 

Thev define a, as the dependence of the absorption coefficient « on the angle of in- 
cidence and it is given by o, = I—|(Z/pc) cos 0—-1/(Z/pc) cos 0+1]?, where Z and pc are 
the acoustic impedance of the material and the medtum, respectively, and @ the angle 
of incidence. Furthermore the average absorption coefficient is given by 
o= eas cos @dé. 

This is a sunulation of a two-dimensional enclosure, although its results permit some 
considerations about the physical problems connected to the reverberation phenomenon. 
It supports the previous authors, who states that the reverberation tume depends on the 
distribution of absorption materials. Even in perfectly diffuse fields, the reverberation 
lume varies as much as 20% with respect to the Sabine’s formula. 

Conclusion is that the angular-dependent sound absorption coeflicient (a, increases 
as the incidence angle increases) has influence on the presence of specular reflections, 
What points to the believing that the sound absorption coefficient possesses the same 


property. 


FE. OTHER LITERATURE 

All the other literature reviewed presented small or null interest for the problem. 
Most of them refer to work in air only and the approximation method do not consider 
normal mode theory and their application to experimental techniques is not convenient. 

From this literature the major points to take is the fact of the reverberation trme of 
cavities With Walls covered by absorption material be related to the acoustical impedance 
of the walls. All the other physical properties of interest, reflection coeflicient and ab- 
sorption can be calculated from it. The article on the measurements of the reverberation 
time in air and water gives some feeling for the numerical results that one may expect 


from experiment. 


All the mathematics of the theories was mnentionally onutted from this shoit :eview 
Because 1 15 of no interest for tls work, all that we were looking for was tlic concepts 


and the mechanics of the problem. 


Hl. THEORETICAL MODEL 


A model for calculating the complex impedance of a low-absorption wall from 
measured reverberation times will be developed. It is assumed that the reader is familiar 
With normal modes in rectangular cavities; otherwise, section 9.7 of KFCS [Ref. 2: 
pp.214-216] provides the necessary background to follow this developinent. All the sub- 
scripts (/, 21, 7) are dropped from the equations for simplicity, but their constant presence 


must not be forgotten. 


The geometry of the problein is presented in Fig. 1. 





Figure 1. Coordinate system. 


A. REVERBERATION IN CAVITIES 

To measure the reverberation properties in a cavity, it needs to be completely filled 
with sound, 1.c., the sound field must be diffuse; the more normal modes excited, the 
better, so that a standing wave pattern is not identifiable. 

The number of normal modes WN that are excited at a specific frequency f, is given 
bv Knudsen [Ref. 3: p.136]: 


poe de Xs | 
vear( =) (1 


Here Tis the volume of the cavity and ¢ the speed of sound in the mediurn. 
The equation governing the sound decav in this cavity with the sound turned off at 


t=O 1s 


es 


p=p(Qe "* Pa 


where t, is the decay factor. 
The reverberation time (7) is defined as the time required for the sound pressure 
level to drop bv 60 decibels (dB), so 


| Pp 
100 — 0 3 
Vlog m(0) nt [3] 
or 
i —T/t 
Oleg = Sloe 2 ‘= —60 [4] 
PO) 


which gives us, after some algebraic manipulation 








i= loa en 
B. DAMPED NORMAL MODES 
Mine Toss Wave Cduatiol 1s 
ae 
Vp= : ao 
c cl 


where p is the complex acoustic pressure and c the complex speed of sound in the me- 
dium, which turn out to be real since we will ignore losses on the medium (relaxation 
time negligible). 


mic Euler equation is 


— 


Cu 
—-Vp= Pp, a [7] 


where p, is the density of the medium and wu the particle velocity. 


~J 


AS @ solution to (6) try 


p= Pcosk,x+,) cos(k,y+o,) cos(kz+ oe (8) 


() 
Where k =—— and the components of the complex wave number are: 


lie = Sage [9] 
ky = Kaales [10] 
ki = k Aye, [11] 


where k,k,k, are the propagation terms and o, o, o, the absorption terms resulting from 


the wall expressed as a bulk absorptron. The complex frequency is written 
oO = wtjf [12] 


where f 1s the temporal absorption coefficient. 


Substituting (§) tnto (6) gives 





= [13] 
or equating real parts 


Aythy +h; —(a;tay+02) = (-Z [14] 


-_, 

|S 
ace 

ie 
Yas 
3s 
a fag 

th 


and imaginary parts 


op 
waa 


6 


Ok tok tok, = 


[15] 


The particle velocity can be found by substituting the value of p , from (8) into (7). 


Its components are: 


ky. . 

ie =i 7a sinfk. x+@,) cos(k,v+,) cos(k.zt+h Jem [16] 
ks, 

ee) as pp cos(k, xt+¢,) sin(kyy+dy) cos(k,z+¢, ae" [17] 


CO 


» Kk, , 7 | 
les cos(kvtd,) cosh v+d,) sin(kz+,)e [18] 


fio : 





= J 


ee gaa VP LICATION 
pcsumie tiacithere is mo-bulkeapserptom-and the wall at ~— J, isithe*only loss. and 
non-rigid wall. This implics @,=0 and ¢, = «,=90. Equation (15) simplifies to 


gp [19] 


5 
a 


C 





Assume that the physical property which determines the absorption characteristics 
of the wall is the normal specific acoustic impedance, (there are no shear waves on the 


walls) and let this be defined as 
Zwall > (Petree [20] 


where y, and y, are normal specific acoustic resistance and reactance coeffierent. 
The normal specific acoustic Impedance at the wall (sav at x = L,) is by definition, 
the ratio of the acoustic pressure (p) to the normal component of the particle velocity 


vector (14), both evaluated at the wall. 


r 
i 
ha 


[21] 


=] 
er * 


At the wall, the wave and the wall rmpedances must match. This gives us the desired 


boundary condition 





x 


—_~ a em: 
Le x= fl “wall aa pos 


or, substituting p (8), 1 (16) and z,, (20) 


Poe INL) = (vt Be es) 
nal oe 


or rearranging. 


() £p 
er er ar 


ee (24 
Vote i Ca 


tan sie) teal sa 


we have the fundamental relation for the physical quantities important to this problem. 


The trigonometric identities 
sin(a+b) = sin @ cos b+ cos asin bh (2s 
cos(a+b) = cos acos b— sin asin b [26] 
and defining’? = x-L, and 0= aL, , casts equation (22) tnto the form 


sin @,, cos(jd)+ cos 6, sin(jo) 


tanlG ot JO) saeeigEcet (a acm) 


[27] 


If 61 (1.e., low absorption) use of (19), and manipulation gives 


: () ; p Oy i og) 0 ctn a 
Y Ol etey Ee i uf eae )| ~— ctn @,-j6 28) 
Working the right hand side separately, multiplying it by (tan @,/ tan 0,) and by the 


complex conjugate of the denommator 


It+/jdctn @, tan Pileac) — O(l+ tan‘é ) 
3 oh ee) ee [29} 
J-—jotan'é, l+0° tan’@, | Fostan cr 


Wilerewsince o-<) the term (1—o-) ai. 


For the left hand side of (28) multiply by the complex conjugate of the denominator 


: CG) : p Ay 
ee | 1 | 
/ Keele) iis) i w Ke )| 


[30] 
keotvetng) US LO be PLN ke] YS 


Now introduce the incidence angle (measured from the normal) w, 





rad 


= 


cos Wf, = . (31] 
From equation (19) and knowing that w/c =k 
k_ B 
a. = ZS [32] 





eh: Oy os ane ee 
Sa SHCUTINS C1) Ga — i 


x 


o,f ” 

— = —— tan’, [33] 
e * e . (1) f 

So the right end side of (28) becomes (using the relation 4, = a cos y,) 


yD) ( fi 2 : fs 2 ie, 
Ses : ne?s{ ™ ban in) vet % tan W)| [34] 
cos wlyy.+7,) | 


and equation (29) becomes 


w gee | W) fi? 
a or eELY Ln. 7 tan Ws) bey — SS a Ltn wm tan | a 
COs Wily cen) COs Wh (yy ) 


x 
tan @(1-8”) — OCi+ tan’0.) 


: j a 
14+-5° tan’@, 1+5° tan’, 


[35] 


i. SOLUTION (general) 
Equating real and unaginary parts of (35), considering 0<1, and (30), we get two 


equations for two unknowns y, and 4,: 


p oe) tan @, 7 
ra son (1. engl) i) oer gee a [36] 
cos Wily ty) Pee Orta a 

and 
ey 2 

| fp ) o(14+ tan’é@_) : 
ee es (15-7, ae Wn) a a LS 7] 
cos WA .+y75) 1+0° tan’, 


Equations (36) and (37) can be more useful if manipulated. First divide one by 


the other and cross multiply 


p 
, O( 1+ tan’é,)— are cau tanw; 


“xX 


= 5) coe [38] 
Ix p 2 2 
tan 0,40 zy tan’ {1+ tan’d,) 


The second equation can be obtained by the subtraction between (36) and (37): 


(146° can*@ 1, (1 i tan’y)—7,(+ is tan’ H) | = 


cos W(y,4+7%) [ tan 0,-S(1+ tan’6,)] 


[39] 


Now define a couple of parameters in order to simplify the expressions (38) and 
(39) for further work (note that I/ tan 6, = ctn @,) 


f 


eee 2 
§ = iS tan Vi; 


[40] 
g = 0(ctn @,+ tan 6,) 


As one may notice, “q” 1s mamly a function of the rigidness of the wall through 
the @, and “s” a function of the incidence angle y, . 
Remember that, so far, the only assumption made 1s “low absorption” of the 


wall (6<1). Let us see what that sav about f/m. Start with equations (19) and (31) 
pee ~ eo [41] 


which, dividing both sides by wm, multiplying the right end side by L,/L,. and remember- 


mg the definiuon of b=a,L, , becomes 


O 
Rie 





COS YW; [42] 


pi 
Ww 
SINCE GOS Ve 1k aldo. = 0 el 1S Clear Wilat pw) 

From this result we may mduce that, if we avoid large vatucs of w,, than s<l , 
1.e., the angle of incrdence can not be close to grazing. 

Fig. 2 shows that it is also acceptable to consider g as being of order 6 if we 
exclude extreme values of ( ctn @,+ tan @,) , 1.e., 8, can not be close to 0 (quasi-rigid wall 


boundary condition) or to z (quasi-pressure release wall boundary condition). 


Pilot of tan(x)+ctn(x 
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Figure 2. Value of “q”. 


It can also be shown that, with some exceptions, s< gq, so that 


ae tan’ YW, < 6(ctn 6,+ tan 6,) 


Using (42) and substituting k by w/c gives 





p OC 
a Stays cos Wy; 
which when substituted in (43) and dividing both sides by 6 gives 


c 


wL, 





<(ctn 6,+ tan @,) 


Assuming that kL,>1 which excludes tow frequencies (already excluded if we wish lots 
of normal modes excited in order to get a diffuse sound field) and_ since 
(ctn 0,4 tan 0,) > 2 (see Figure 2), the relation is true if we avoid grazing incidence (al- 
ready excluded). 


The assumptions and restrictions, so far, are: 


- low absorption at the wall (6<1) 
- exclude close to grazing mcidence (79a) 
- avoid rigid or pressure release wall, (0,2 0 and 0, 4 2/2) 


eS 1G le 


Note that the validity of the last relation can be extended when the dimensions 
of the cavity or the frequency are increased. 
2. RESULTS 
In this sub-section, expressions will be found for the normal resistance and 
reactance of the wall, using the normal incidence case as a check point. Insertion of the 


parameters s and g into (38) and (39) yields 








ve g= 7s 
Wx  i+sq ine] 
nx(I—s) — y,(1-+s) = (y¢+n%) cos wy, tan 6,(1 — q) [47] 


Since ox<! anditang. 1s net laree, the tem (1-07 tat. 

Applying the approxmnations s < g<lI, which implies (q—s)<l and sgl, it can 
be seen that the denominator of (46) goes to I and that (1+s)=(l—s)=(l—g)=1 The ratio 
/#, 18 much less than |. It also tells us that both the resistance and the reactance have 
the samme sign. 


Equations (46) and (47) become: 


Fe ages [48] 
Nx—Voxlystns) cos w; tan @, [49] 


Taking the value for y, from (48) and substituting on (49) yields: 
nxl1—(q—s)] = [(q~-s)' +1] cos ; tan 8, [50] 


Dividing both sides by 7, and noticing that (g—-s)<1 , we have an expression for ,: 


| 
cos w; tan @, 


[51] 


ae 


Now substitute 7, fron) (51) into (47) and get y, 


G—S$ 
eC Om atc 


or, plugging back the vatues for the parameters s and q (39): 


5 tan’, 
—-— ( ctn’0+ 1)— aks wen [S53] 


y = 
tor CCONNT, OY Cosy tan 0. 


The angle of incidence varices from 0 to z/2, so cos y, 1s always positive, @, 1s a 
function of the mode number and the “rigidness” of the walls, so 1t can take any value, 
1.e., tan @, can be either positive or negative. Under these arguments, looking to (51) and 
(52), and remembering that (g—s) > 0 it 15 easy to verifV that if tan @, >0, then both y, 
and #7, are positive and they are negative if tan @, < 0. This is consistent with the expec- 
tations Stated earlier in this section. 


Applving the relation 4, = w/ecos w, to @, vields 





Een . 
== COS Yi; [Ss] 
ondecs can berelated to ff by (31) 
ead : 
Oo. = ao 
oe OS FCG ee 


The temporal absorption coefficient f can be obtained from the reverberation 


time 7 by direct comparison between equations (2) and (9) (on its maximum): 
ee —pt 2 ayo? 5 
Pp = pplOje © = pp(DJe °* [56] 


which leads to 





p= (57] 


and with t, and 7 related by equation (5), we have 


= 
\| 
———s 
| ge) ae 
“len 
[it 
tN 
cif 
LJ) 





— 
sr 


So al] the variables in (83) and (51) are known or can be found from measure- 
ments, Which permits the calculation of y, and 9, of the wall. 
52 (CHEGES 
a. Alain assumption 
First we can check the assumption 0<I, to do so we may solve (51) and (53) 


for 0 


ye ie 


(ctn @,+ tan @,)— 


ee - [59] 


wL,tany; 
The denominator is always greater than 0 (except for grazing incidence) , and it is larger 


a¥, 
/ 


than 1 because we assumed it when checking for s < q, since ae <1 our initial assump- 
tion is verified. 
b. Normal incidence 
We can start with (38) and (39), considering again that (1+0? tan?@é,)~1 and 
w,= 0. Without any other assuinpuions that 6<1 and @,%(7+1)z/2 or nz, equation (38) 


becomes 


vy (14+ tan’@,) 








Wx tan, LOM) 
SUBSILULION Of ys mi (5 oyvicids 
nd 1~q) = 93[ 1+(1—q)"] tan 6, [61] 
dividing both sides by 4, and noting that (l—qg)~I we get 
eee [62] 
batho. 
which, substituting in (60) gives 
q 2 
a 6, = 0(ctn’é,+1) [63] 


This results are the same as (52) and (53) if we set yy, =0. 
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¢. Power reflection coefficient 
Nespower tenection eoellicicnt (i) 1s given ny KICS jer 21 139] by 


equation (6.45) 


PR (y, COS ae Pn, cos’ Ws [64] 
° ' 2 i 2 
‘ (7, cos W+1)"+7, cos yy; 


We want it to be close to one (low absorption), so, dividing evervthing by ( cos?y,) 
Aree | ere : 
(y,— sec ),)}) +a, =(yet sec W)) +7), [65] 


This tells us that 7, can take any value, and 


y—2y, sec yt sec w= t2y, sec Wj+ sec? Wi [66] 
O1 
Vx SCC Hi=y, SCC Y; [67] 
which becomes 
Y~€ COS [08] 


l'rom this relation we find that y, despite of being much smaller than 7, has 


to be small. Recalling equation (46) one may see that 


qs 





) — 


oS T+sq HK COS W; [69] 


which, after some manipulation, and taking into account the approximations earlier 
made, tells us that 7, has to be also much smaller than (cos ¥,). So we expect, for low 


attenuation at the walls, low resistance and low reactance. 


LV CONE USO IS 


Using normal mode theory in cavities it 1s possible to calculate the normal acoustic 
resistance and reactance as functions of the reverberation time, the frequency of the 
acoustic field, the speed of sound in the medium, the dimensions of the cavity, and the 
incidence angle. 

To measure the reverberation time, a diffuse field is required, t.e., a fairly uniform 
pressure is desired throughout the cavity. Tle reverberation time can be related to the 
temporal absorption coeflicient and it is expected to increase with the increasing di- 
mensions of the cavity, the decreasing absorbtion of the walls, and the decreasing pc of 
the medium. 

Under the following hmitations: 

- low absorptive wall 
- avoid rigid or pressure release walls 


- exclude angles of incidence close to grazing 


the results obtained as a function of deternunable parameters are: from (53), the normal 


specific acoustic resistance is 








2) 
p a eek a) pi tan y; a 
eee PCW ea COC ict al ah eG [70] 
ccos W; io 
COSY ttm a 
C COS W; 


from (51), the norma! specific acoustic reactance Is 


[71] 


eS 





ieec) 
cos W;tan{ —2— cos y; 


It 1s noticeable that, under the mentioned limitations, the resistive part has to be 
small and much smaller than the reactive part. 

The exclusion of the “close to grazing” angles of incidence is assurned as not being 
very important because, for a diffuse field, there are many excited normal modes, so the 


modes excluded by this limitation are relatively few. 
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Despite the limitations on this model, it is a gencralizauon of the theory presented 


ii the textbooks | Reis 3. Gi. 


to 
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